Semiclassical theory oi h/e Aharonov-Bohm oscillation 
for doubly connected ballistic cavities 
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In Aharonov-Bohm (AB) cavities forming doubly connected baUistic structures, h/e AB oscil- 
lations that result from the interference among the complicated trapped paths in the cavity can 
be described by the framework of the semiclassical theory. We derive formulas of the correlation 
function C{/S.(f)) of the nonaveraged magnetoconductance for chaotic and regular AB cavities. The 
different higher harmonics behaviors for C(A(/)) are related to the differing distribution of classical 
dwelling times. The AB oscillation in ballistic regimes provides an experimental probe of quantum 
signatures of classical chaotic and regular dynamics. 

PACS numbers; 05.45. +b, 03.65.Sq, 72.20.My, 73.20.Fz 

Electron transport through quantum cavities is an exceedingly rich experimental system, bearing the quantum 
signature of chaos. [Q On the theoretical form, powerful techniques based on semiclassical approaches have produced 
specific predictions testable by experiments. 1^-^ An interesting result that has emerged concerns the magnetotrans- 
port of doubly connected ballistic cavities, i.e., Aharonov-Bohm (AB) cavities ]7|-[To|] (see Fig. 1). We have calculated 
the average conductance for these systems and showed that the self-averaging effect causes the h/2e Altshuler-Aronov- 
Spivak (AAS) oscillation, which is ascribed to interference between time-reversed coherent back-scattering classical 
trajectories. Moreover we have showed that the AAS oscillation in these systems becomes an experimental probe of 
the quantum chaos. Another interesting phenomenon in these systems is the h/e AB oscillation for nonaveraged 
conductance. The result of numerical calculations |q] indicated that the period of the energy averaged conductance, 
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changed from h/2e to h/e, when the range of energy average Ai? is decreased. However, little is known about the effect 
of chaos on the h/e AB oscillation in AB cavities. In this paper, we shall calculate the correlation function C(A(/)) 
of the nonaveraged conductance by using the semiclassical theory and show that C(A0) is qualitatively different 
between chaotic and regular AB cavities. 

In the following, we shall derive C(A0) separately for chaotic and regular AB cavities in which uniform normal 
magnetic field B (AB flux) penetrates only through the hollow. 
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FIG. 1. An example of the pair of four classical paths which contribute to the correlation function for the Aharonov-Bohm 
cavity. 
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The transmission amplitude from a mode m on the left to a mode n on the right for electrons at the Fermi energy is 
given by 



tn.rn = -ifi^VnVm dy dy' il)*^{y')'4>m{y)G{y' , y , Ep) , 
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where Vm{vn) and tpmi'^Pn) are the longitudinal velocity and transverse wave function for the mode m (n) at a pair 
of lead wires attached to the billiards. In eq. (|^), G is the retarded Green's function. In order to carry out the 
semiclassical approximation, we replace G by the semiclassical Green function, 
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where Ss is the action integral along a classical path s, the pre-exponential factor is 
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with 9 and 9' the incoming and outgoing angles, respectively, and /i is the Maslov index. Substitute eq. into eq. 
and carrying out the double integrals by the saddle-point approximation, we obtain 
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where W is the width of the hard-wall leads and m = ±m. In eq. (||), Ss{n,m; E) = Ss{y'o,yQ', E) + hTr{myo — 
ny'„)/W, Ds = {'m,VFCOs9')-^\{dyld9')e\ and /I, = ^is + H {-{89 / dyYy) + H {-{89' /dy')g) , respectively, where 
9 — s\n^^{nTT/kW) and H is the Heaviside step function. 

Transmission coefficients between modes are obtained by taking the absolute square of transmission amplitudes, 
. For leads of width W that support Nm = Int[/cM^/7r] modes, the total transmitted intensity summed 



over m and n is 



T{k) 



1 TT 

2kW 



Nm 



(6) 



where s and u label the classical trajectories. In eq. (^, Lg = Ss/kfi , Vs,u = (A« ~ Ms) /2 , and Ag = {fik/W) Dg- 
The fluctuations of the conductance g = {e^ / 'itK)T {k) are defined by their deviation from the classical value; in the 
absence of any symmetries. 
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In this equation gd = {e^ /ttK)Tci, where Td is the classical total transmitted intensity. In order to characterize the 
h/e AB oscillation, we define the correlation function of the oscillation in magnetic field B by the average over J5, 

C{^B)^{5g{B)5g{B + ^B))^. (8) 

With use of the ergodic hypothesis, B averaging can be replaced by the k averaging, i.e., 

C{^B)^{5g{k,B)5g{k,B + ^B))^. (9) 

Within the diagonal approximation the correlation function of transmission coefficients between the modes is 
given by 



Cd{^B) = 
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where J^n m ^'^n,m. = T{k) — Tel- The semiclassical expression for the transmission amplitudes, eq. (||), yields 
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X < exp 
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where 9 — ±9. As for AAS oscillation, the diagonal approximation yields an expression with k dependence only 
in the exponent. With use of Ss{B) — hkLg + e A • dr, we get 
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Here 2tt J^^^^ A-dr — iJG^^^) and 2tt J^^.^^ A' - dr = {B + AB)Q^yy The finite k average implies that only contribution 
is expected for 
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exactly. Because of the definition of C_d in cq. all four paths satisfy the same boundary conditions for angles, 
and hence they are all chosen from the same discrete set of paths. In the absence of symmetry, the only contribution 
is w = s and t = u. The terms with s = u and t = v are excluded because they represent the average values that 
must be removed from the correlation functions. In Fig. 1 we show the typical set of trajectories that contribute to 
the correlation function. This process is analogous to the two diffuson propagators in a diffusive regime jl^]. Since 
magnetic flux penetrates only through the hollow, the exponent in eq. ( |l2|) becomes 
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where ± corresponds to the clockwise (counterclockwise) rotation to the center disk for path u. In eq. ([Tj) Wa is the 
winding number of classical path s. Therefore we obtain 
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In order to evaluate sum over s and integrations on 0{9'), we shall reorder the trajectories according to the increasing 
dwelling time Tg . Therefore we find for the diagonal part of the semiclassical correlation function for chaotic systems 
as 
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where 5 = y^zTVy/a. In deriving eq. ( |l6|) we have used the exponential dwelling time distribution, N{T) ~ exp(— 7r), 
and the Gaussian winding number distribution for fixed T, [p^ i.e.. 
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where Tq and a are the system-dependent constants corresponding to the dwelling time for the shortest classical 
winding trajectory and the variance of the distribution of w, respectively. By using the extended semiclassical theory, 
pM we can take account of the off-diagonal part and the influence of the small-angle diffraction as 
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for the case which the widths of the lead wires are equal. Then we obtain the full correlation function for chaotic AB 
cavities, 
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The periodic function C(A0) has the minimum value 
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at A(f>min = niT, where n = 1,3,5, •• •. Therefore, C{A(p) oscihates with the period (/)o, i-e., AB oscillation. From 
the above results, we can conclude that it is possible to predict quantitatively C{A(f>) of the chaotic AB cavities from 
a knowledge of the chaotic classical scatterings dynamics. Note for consistency that the field scale of fluctuations 
is twice that of AAS oscillation |Q because the relevant phase involves the difference between two winding numbers 
whereas AAS oscillation involves the sum. Surprisingly the semiclassical formula, eq. (p^), is quite similar to Isawa 
et al.'s results for the disordered quasi-one dimensional AB ring: [T^ 
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In this equation is the phase coherence length, R is the radius of the ring and Q = smh{2TTR/ L^)+suih^ (ttR/ L^)/2, 
respectively. 

The periodic function C{A(f)) is large and positive for very small A(/), and has the limiting value 



(22) 



This result is consistent with the result of random matrix theory for the circular orthogonal ensemble. ||T^,^ In the 
case of weak A(j), eq. (|l^) is rewritten asymptotically as 
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Therefore C(A(^) decreases quadratically with increasing Acf) near A(f) = 0. The quadratic behavior of C(A(/)) is 
similar to that for ordinal chaotic cavity, e.g., stadium, at near AB = 0. |^,|| 

On the other hand, for the regular cases, we use N(T) ~ [ppc|] in eq. ( p^ Assuming as well the Gaussian 
distribution of P{w;T), we get 



C{A4>) = C(0) COS ( 27r^^ { 
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(24) 



where F is the hyper-geometric function of confluent type. As in AAS oscillation, parameters (5 and a characterizing 
the classical dynamics determine the behavior of C(A(/)). 

Next we shall see the difference of C(A(/)) for chaotic and regular AB cavities in detail. In the chaotic AB cavity, a 
main contribution to the AB oscillation comes from the n—\ component. Figure 2 shows (a) aspect ratio [a — R/W) 
and (b) the degree of opening to the lead wires (77 = L/W) dependence of C{A(p) for the open chaotic AB cavity 
(Sinai billiard |^l|), where R is the radius of the center circle and L is the linear dimension of the outer square. 
The classical parameter S is calculated by using geometric dimensions of the cavity. ||^,^ In the case of small a or 77, 
classical trajectories are able to wind around the center disk many times and the higher harmonics can contribute to 
AB oscillation. Therefore one can see from Fig. 2 that the minimum value C (A(f>min) / C (0) slightly increase from -1 
as (7 or /X becomes small. 
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(a) (b) 




FIG. 2. Semiclassical correlation function C(A(j!)) of the conductance as a function of for the chaotic AB (Sinai) cavity: 
(a) for various aspect ratio a = R/L {rj = 0.2) ; (b) for various degrees of opening rj = W/L {a = 0.3). 

On the other hand, for regular cases, the amphtude of the AB oscillation decays algebraically, i.e., F ^ n^"^^^^ for 
large n. This behavior is caused by the power law dwelling time distribution, i.e., N{T) ~ T^^ . Thus, in contrast to 
the chaotic cases, we can expect that the considerably higher harmonics contribution causes a noticeable deviation 
from the cosine function for C(A0). Therefore, between the difference C(A(/)) of these ballistic AB cavities can be 
attributed to the difference of chaotic and regular classical scattering dynamics. 

In summary, we have investigated magnetotransport in single ballistic cavities whose structures form AB geometry 
by use of semiclassical methods with a particular emphasis on the derivation of the semiclassical formulas. The 
existence of the AB oscillation of nonaveraged magnetoconductance is predicted for single chaotic and regular AB 
cavities. Furthermore, we find that the difference between classical dynamics leads to qualitatively different behaviors 
for the correlation function. The AB oscillation in the ballistic regime will provide a new experimental testing ground 
for exploring quantum chaos. 

We would like to acknowledge K. Nakamura and Y. Takane for valuable discussions and comments. 
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